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Abstract 

A precise definition of an adiabaticity parameter v of a time-dependent Hamiltonian 
is proposed. A variation of the time-dependent perturbation theory is presented which 
yields a series expansion of the evolution operator U (r) = J2e U ( r ) with (r) 
being at least of the order z/. In particular U^°\t) corresponds to the adiabatic 
approximation and yields Berry's adiabatic phase. It is shown that this series expansion 
has nothing to do with the 1/r-expansion of U(r). It is also shown that the non- 
adiabatic part of the evolution operator is generated by a transformed Hamiltonian 
which is off-diagonal in the eigenbasis of the initial Hamiltonian. Some related issues 
concerning the geometric phase are also discussed. 
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1 Introduction 



Since the introduction of the adiabatic geometrical phase by Berry [jl] , the study of the cyclic 
evolution of quantum states of non- conservative quantum systems (explicitly time-dependent 
Hamiltonians) has attracted much attention. By definition a pure cyclic state |at; t)(jv; t\ is 
an eigenstate of the evolution operator U(t), i.e., 

U(t) =: £ e^ (r) k; r) r| , ^(r) G [0, 2tt) . (1) 

Therefore the quantity of main importance in the study of cyclic states and the accompanying 
phases /^(r) is the evolution operator U(r). This operator is defined by 

|V>(r)> =: I7(r) |V>(0)> , (2) 

where \ip{r)) is the solution of the Schrodinger equation. Alternatively, one can define U(r) 
as the solution of 

±V(r) = -\H [T )U ( r) (3) 
(7(0) = 1, 

where if = if (t) stands for the Hamiltonian and r G [0, oo) parameterizes the time. The 
solution of Eq. (§) can be implicitly expressed in the form: 

(/(r):=Te-^o H ^, (4) 

where T denotes the time-ordering operator 0. 

The purpose of this article is two fold. First, a precise definition of an adiabaticity 
parameter v will be given. This is a parameter which quantifies the rapidity of the time- 
dependence of the Hamiltonian. Next, a series expansion of the evolution operator U{r) 
will be proposed whose £-th term is at least of the order iA In particular the adiabatic 
approximation: 

u(t) «xy Q " (T) l n ; r )( n ;°l , (5) 

n 

is recovered as the first term. In Eq. ([|), \n;t) are the instantaneous eigenstate vectors of 
the Hamiltonian, i.e., 

H(t)\n;t) = E n (t) \n;t) , (6) 
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and 

a n (t) := 5 n (t) + ln (t) , 5 n (t) := -~ J* E n (t')dt' , ^y n (t) := i j\n;t'\ — \n;t') dt' . (7) 

If the Hamiltonian is T-periodical, i.e., H(T + 1) — H(t), the a n (T), S n (T), and ~f n (T) are 
called adiabatic total phase angle, adiabatic dynamical phase angle, and adiabatic geometrical 
phase or Berry phase angle, respectively [1J. 

The approach pursued in this article clarifies the misconception that the adiabatic ap- 
proximation of U(t) is given by the zero-th order term in its 1/r-expansion, i.e., it is only 
valid for r — > oo. It also provides a convenient framework to perturbatively compute the 
non-adiabatic corrections to Berry's adiabatic phase. 

2 The Adiabaticity Parameter 

Following Berry consider a parameter-dependent Hamiltonian: 

H[R] = y £E n [R] \n;R)(n;R\ , (8) 

n 

where R = (R l , ■ ■ ■ , R d ) are real parameters viewed as the coordinates of a parameter space 
M., E n [R] are the real eigenvalues and {|n; R}} is a complete orthonormal set of eigenvectors 
of H[R\. In this case a time-dependent Hamiltonian H(t) corresponds to a curve C : [0, r] — > 
M in M: 

H(t) := H[C(t)] = H[Rc(t)) , with C(t) =: ■ ■ • , R&(tj) =: R c (t) . 

Furthermore, suppose that the spectrum of H[R] for all R G M. is discrete, the energy 
eigenvalues E n [R] are non-degenerate, and there is no level- crossing, i.e., for every t G [0, r], 
E n [R c {t)} < E n+1 [R c (t)]. 

In order to define a dimensionless adiabaticity parameter, first one defines a characteristic 
frequency: 

w c (ti,t 2 ) := Sup{|4wi(*)| ■ n^m = 0,l,--- , te[T 1 ,r 2 ]}; [n, r 2 ] C [0, r] , (9) 
with "Sup" abbreviating supremum (least upper bound) and 

A mn {t) := (m;t\^-\n;t) (10) 
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m;t\ f t H(t) \n;t) 



for m 7^ n . (11) 



E n {t) — E m {t) 

Then the desired adiabaticity parameter is defined according to 

" = -^KF 1 > (12) 

where AE is a convenient energy scale. For example one can take AE to be the first 
transition energy of the initial Hamiltonian, i.e., AE := Ei(0) — Eq(0). 

The main motivation for this definition is the fact that according to Eqs. (|9]), (|IID , and 
(|I2|), uj c (0,t) and therefore v involve time derivatives of the Hamiltonian H(t). The use of 
u; c (0, r) in the definition of the adiabaticity parameter will be self-evident once one examines 
the evolution operator. The role of AE is to provide a convenient energy (fox frequency) 
scale. 

It must also be emphasized that by definition, v is a "global" quantity which characterizes 
the time-dependence of the Hamiltonian. In particular, it is neither equal nor proportional 
to 1/r. Furthermore, note that although \n;t) are only determined up to arbitrary Rc{t)- 
dependent phase factors (gauge transformations along C, @), |A mn (t)| with m ^ n and 
therefore a; c (ri,r 2 ) and v are independent of the choice of such phases (they are gauge- 
invariant quantities). In other words, v is well-defined. 



3 Computation of U (r 

Consider the definition of the time-ordered product in (||) 

U(r) 



lim 



1 - -H(t N )e 



1 - -H(t k )e 



1 - j-H(t )e 



(13) 



where e := t/N and t k '■= ke, for k = 0, 1, • • • N. Using the orthonormality and completeness 
of the energy eigenvectors |n; t), one can compute: 



U(t) 



no,---,n N 



lim \\nN;t N )(n N ;t N \ 



1 - -H(t N )e 
n 



l^fcj tk){Uk\ tk\ 



1 - -H(t k )e 
n 



n>N-i] tN-x) (njv-i; tjv— 1| 

n k-l] tk-l)( n k-l] tk-l\ 
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|ni;ti)(ni;ti 



1 - j-H(t )e 



\n ;t ){n ;t \ 

N 



/ N N 

lim \e-^=° En ^ )e U(n k ;t k \n k ^;t k 

N — >oo \ 

no,—,n N \ k =l 



-1) \n N ;t N )(n ;t \ . 



Introducing 



(n N ; T\U(r)\n Q ] 0) 



lim 



N 



Y[(n k ; t k \n k -i] t k -i) 

k=l 



one then has: 



U{t) = K mn {r) \m;r)(n;0\. 

mn 

The computation of the terms in the product in Eq. (0), is straightforward: 



(14) 



(15) 



(n k ; t k \n k -i] t k -i) 



+ 0(e 2 ) 



t—tk-i 



where 



(16) 

A mn (t) . (5mn l)-^mn(^) • (^v) 

Substituting ([16]) in (0), one obtains 2^ terms which can be arranged in the order of 
the appearance of different powers of e. In this way one finds only N terms of order e° = l.Q 
These will be denoted by i^^(r): 



iV 



iV— »oo 



T . 



1=0 



Performing the algebra one finds: 



K® (r) 



E 


E 


e 






£ =0 


^noni 













il+l 



-i(^)^ni n/ + i ' ' ' ^njy.iniv 



+ 0{e) , (19) 



1 Note that a sum of iV^ terms of order e e is of order e° = 1. 
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where £ — 0, 1, • • • , N and t k = ke/N = 0, e, • • • , r. In particular, 

K^lnAr) = T nN (r)5 nNn0 + O(e), 



(20) 



^no(^) = -r * H 'n N n (ti) + 0(e) , 

" +-—4-. 



^4^ 



ti=0 



+ 0(e) 



(21) 



where 



J- n^fc J — e J . 

<„(t fc ) := ihr* m (t k )A mn (t k )r n (t k ) , 



(22) 
(23) 



and r^(tfe) = l/r m (tfc) is the complex conjugate of r m (^). 

Furthermore, for every £ > 1, one can express K$ n (r) in terms of K$ n (r) and H' mn {ti), 
namely 



ifW (r) = — V V tf(°> (r) 

m m e \ J \ fj I ' ' m mo \ I 

V ' t il <-<ti e =Q mo-mn 



IT ^m„-im (^a/ 
Lo=l 



+ O(e) , 



•t . \ it / mo-m^_i 



n H ma _ ima (ti a ) 



.a=l 



+ 0{e) . (24) 



Therefore if one defines: 
C/W(r) 
lf'(f) 



:= lim £tfffi(T)|m;r)<n;0|, 

mn 

■= i im J2 H 'mn(t)\m;0)(n;0\, 

N^oo 

mn 

= ih^e- i[am{t) - an{t)] A mn (t)\m;0)(n;0\ , 

mn 

= J2 e -<[«™(*)-«"(0] A„ m (t) |m; 0)<n; 0| , 



(25) 



then 



U(t) =f^U^(r) =C/(°)(r) 



£=0 



(26) 



(27) 
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where by taking iV — > oo the sums of the form J2t=t e f(ti) have been promoted to the 
integrals // dtf(t). In particular one has: 

[/(°)( T ) = J2e ian(r) \n-T)(n;0\ , (28) 



C/W( r ) = £/(o)(, 



— i 
ft Vo 



(it if'(i) 



(29) 



By construction U^(r) yields U(r) in the adiabatic approximation @. Therefore, adia- 
batic approximation is valid only if one can neglect /J" dtH'(t) in ( 27f) . It is not difficult to 
observe that this condition is fulfilled if the adiabaticity parameter v as defined by Eq. ([T2] ) 
is negligible. In fact, for every t\ and ti satisfying < t\ < £ 2 < r, one has: 

|(m;0| f 2 dtH'(t)\n;0)\< dt\(m;0\H'(t)\n;Q)\ <% dt\A mn (t)\ < At AE v , (30) 
Jt± J ti J ti 

where At := r 2 — T\ with [r 1; r 2 ] C [0, r] is the time interval over which u c is non- vanishing, 
i.e., T\ and r 2 and therefore At are defined by the condition: 

uj c (t[, t' 2 ) = if and only if r 2 < T\ or r{ > r 2 . (31) 

Note that usually At is a finite time interval whereas r may be infinite (arbitrarily large). 
A case for which At = r is when the Hamiltonian depends periodically on time. In this case 
however, the physically interesting features of the evolution is given by At = t <T, where 
T is the period of the Hamiltonian. Particularly interesting is the case r = T. 

In general, At may become arbitrarily large in which case ( |30"D is a trivial statement. 
However, even in this case v plays a most important role. In fact, one can argue that in 
general the following relation holds: 

AT-l 

U{t)=Y j U^\t)+0{u n ). (32) 

In particular the adiabatic approximation is valid if and only if v <C 1. In order to es- 
tablish Eq. (|32D , it is sufficient^ to show that U^'(t) and therefore / dt H'(t) is at least of 
order v. This is however self-evident since one knows that if v — 0, then U(r) = U^(t). 
Consequently, in a //-expansion of U(t), (r) is necessarily of order v or higher. 

2 This is because U™'(t) involves £ copies of H'{t) in its integrand. 
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An alternative way of expressing Eq. (|3"2|) is to define U^(t) := U"'{t)/v and H'{t) 
H'{t)/v and write U(t) in the form: 



e/(t) = jr# w (7-) = 



(33) 



This equation must not however be viewed as a true power series expansion of U(t) in z/, 
for U^'(t) may in general depend on v. 

A simple consequence of Eq. (|33|) is the fact that the essential ingredient which determines 



U(t) beyond the adiabatic approximation is not r, but (ti, T2) of (31) and ia For physically 
realistic non-periodic Hamiltonians, the latter quantities take finite values whereas the du- 
ration r of the evolution of the system may be arbitrarily large. For a quantum system with 
a periodic Hamiltonian H(t) = H(t + T), the physically interesting case is when r 2 = r < T. 
In this case, v may or may not be determined by T. This is because, in general the pa- 
rameters Rc{t) of the Hamiltonian have different periods. These are necessarily of the form 
T M = T/z^, respectively, where z^ is a positive integer. Clearly cj c (0,t) and consequently 
v will depend on the largest value of z M . Therefore, in general the statement that "the 
adiabatic approximation is valid if the period of the Hamiltonian is large," is false. The 
only case where v depends on T and the preceding statement is valid, is the case where the 
energy eigenstates are time- dependent and either there is effectively one changing parameter 
or T M = T for all fi. Typical examples of these two cases are a spin in a precessing magnetic 
field and a spin in a precessing and nutating magnetic field with equal precession and 
nutation periods, respectively. 

Another useful observation is that in the eigenbasis {\n; 0)} of H(0), H'(t) has no diagonal 
matrix elements, i.e., 

{0 for m = n 

(m;tU.m)]\n;t) i^(t)-ajt)] r , ( 34 ) 

ltl E m (i)-E n {t) e loi m^n, 

where use is made of Eqs. (53) and (fLip. This equation also implies that the matrix elements 
(m; t\U^> (t)\ti\ 0) are directly related to ( f~^r J > i- e -> the main contribution to U^(t) comes 
from the nearest energy levels. 
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4 Relation between H'(t) and H(t) 



Consider Eq. (P^). Since the m = n term is missing in the sum in (PEj), one can express 
H'(t) in the form: 

F'(t) = -iU J2 [(m-t\e- iam{t) \ j t [e iaJt) \n;t)] |m;0)(n;0| . (35) 

This equation together with the definition of U^°\t) leads to 

H'{t) = U m {t)H{t)U {0 \t)-i?iU m {t)j t U {0) {t) . (36) 

The resemblance of this equation to the gauge transformations of the gauge potential in non- 
Abelian gauge theories (connections on principal fiber bundles) is remarkable, [[| . Moreover, 
a simple calculation shows that 

U<M(t)H(t)UV>(t) = 5X(*)|n;°><n;°l. (37) 

n 

H'(t) = -i% J2( m ^\U m (t)j t U m (t)\n;0) \m;0)(n;0\ . (38) 



Therefore, the role of the first term on the right hand side of (|36| ) is to cancel the diagonal 
matrix elements of the second term. 

Eqs. (^) and (|38| ) signify the importance of the adiabatic approximation in the deter- 
mination of the non-adiabatic corrections U^(t) (£ > 0) to U^(t) and consequently the 
exact evolution operator U(r). 

In order to fully appreciate the meaning of Eq. (|36|) , one should recall the effect of a 
general time-dependent unitary transformation of the Hilbert space Ti. Let tl(t) : 7i — >• 7i 
be such a transformation and suppose that the transformed state vectors 

\#(t)):=U{t)W)), (39) 
satisfy the following Schrodinger equation: 

l nj t \m) = H{t)m)). (40) 
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Then requiring \i(^(t)) to satisfy the Schrodinger equation defined by the original Hamiltonian 
H(t), one finds: 

H(t) = U(t) H{t) U\t) - i% U(t) j U\t) . (41) 

Therefore Eq. fl36|) is a particular example of Eq. (gl]) with U{t) = U m (t) and H(t) = H'(t). 
In other words, the non-adiabatic effects are given by a transformed Hamiltonian H'(t), where 
the transformation is performed by the inverse of the adiabatically approximate evolution 
operator, i.e., U^(t). 

Transformation (|39|) may be used to yield an exact solution of the Schrodinger equation, 
if one can find a lA{t) that makes the transformed Schrodinger equation (|40| ) easily solvable. 
The Rabi-Ramsey-Schwinger method is a special case of this approach where one uses 
the symmetry of the problem to find a unitary transformation U{t) which renders H time- 
independent. 3 This is in general possible for cranked Hamiltonians with a fixed cranking 
direction || |7|]. 



5 Time-Dependent Unitary Transformations and the 
Geometric Phase 

There is another significant application of time-dependent unitary transformations of the 
Hilbert space when the Hamiltonian is T-periodic. In this case, one knows from the Floquet 
theory || that the evolution operator is given by U(t) = Z{j) exp[— zrif /h], where Z(r) 
is a T-periodic unitary operator with Z(T) = Z(0) = 1, and if is a time-independent 
Hermitian operator. Then by definition the cyclic states \m;T)(m;T\ at r = T are given by 
the eigenstates of H and the associated total phase angles are Pm(T) = —T{m; T\H\m; T)/h. 
Next consider choosing U{t) of Eq. © to be Z\t) and |^(0)) = |^(0)) = \u;T). Then by 
virtue of Eq. (|4l|) , one has: 



PAT) = -~Mo)W(o)) = r 



A M (T)+T M (T), (42) 



3 For a recent review of the application of this method to a spin system in a precessing magnetic field, see 
Ref. §. 
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A^r) := --/ r//-(r(/)|f/-(/)|r(/-)} . ( J3) 



o 

r^(r) := tJ o T dt(Mt)\j t \Mt)), (44) 

where := Z(£)|^(0)). The phase angles A jV (r) and l~V(r) were first introduced by 

Aharonov and Anandan ||, and called the dynamical and geometrical phase angles.0 They 
generalize their adiabatic counterparts previously discovered by Berry IJ. Clearly, for !/<l, 
one has: \m;T) « |n;0), A^(T) « <J n (T), and IV (T) « 7n (T). 

Next consider the general case. If U(t) is also assumed to leave the original energy 
eigenstates invariant, i.e., H(t) has the same eigenstates as H(t), then one can show that 
U(t) and H(t) must be of the form: 

U(t) = e if{t) U , (45) 
H{t) -> ff(f) = if (*) - ft ^ 1 , (46) 

where / = /(i) is a real- valued function, Wo is a constant unitary transformation, and 1 is 
the identity operator on 7i. Note that transformation (|39| ) with W(t) given by ( |4"5| ) leaves all 
the observables invariant. It shifts the energy eigenvalues: 

E n (t) - K(f) = E n (t) - ft ^ , (47) 

but does not affect the observable transition energies E n (t) — E m (t). Therefore, one can 
identify a physical system with the equivalence class of all the quantum systems which are 
related with transformations of the form ( |39| , [45D . 

A simple consequence of this observation is that unlike the dynamical phase angle (j43|). 
the geometric phase angle ([44]) is a physical quantity [0. This statement however requires 
a clarification. One can easily see that ( [46]) changes the dynamical phase angle A AA (r). It 
does nevertheless leave the difference of two dynamical phase angles and consequently two 
total phase angles invariant. Therefore, the dynamical phase angle also carries physically 
significant information about the evolving system. In fact, a close look at the experimen- 
tal results on the detection of the geometric phase fl(| clearly shows that the measurable 



4 Note that in general the Hamiltonian needs not be periodic and the Floquet theory does not apply. 
However, even in this case one can define \<p//(t)) and Eqs. are valid. See fi 0] for more detail. 
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quantities are related to the differences of total phase angles. This is usually overshadowed 
by the fact that in the best studied system, namely the two- level spin (j = 1/2) system in 
a precessing magnetic field |]J, the conventional choice of the Hamiltonian (|48|) leads to two 
cyclic states whose total phase angles differ by a minus sign. Therefore their difference is 
twice one of them and it appears that the experiments detect a total phase. In reality how- 
ever, the experiments always detect the differences between total phase angles, a quantity 
which is invariant under P6|). This can be easily seen if one uses a complete orthonormal 
set of cyclic states {|at; t)} to compute the expectation values: 

<X(r)> := ^(r)\X(0)mr)} = (m\U\r)X(0)U(r)mO)) , 

= J2(mW;r)(^r\U\r)X(0)U(r)W,T)^;r\m) , 

where X(0) is an observable and in the last equality used is made of Eq. ([[]). 

Incidentally, it is not too difficult to see that by definition, the adiabaticity parameter 
v (|I2"D, the operator H'(t) of Eq. and consequently all the non-adiabatic corrections 

C/W(r) (with I > 1) to U^(r) are invariant under the transformation Therefore all of 

these quantities signify physically measurable effects. 

6 Application to a Magnetic Dipole in a Magnetic Field 

In this section, the utility of the expansion (p7|) is demonstrated for the quantum system 
consisting of a magnetic dipole moment subject to a magnetic field with changing direction. 
The HamiltonianQ of this system is given by: 

H(6, <p) = bR(9, tp)-J= 6(sin 9 cos (pJ 1 + sin 9 sin tpj 2 + cos 9 J 3 ) , (48) 

where b is the Larmor frequency, 9 and (p are the azimuthal and polar angles in spherical 
coordinates, respectively, and J is the angular momentum operator with components J M , 
fi — 1, 2, 3. Clearly the parameter manifold is the two-dimensional sphere S* 2 and the time- 
dependence of the Hamiltonian is described by a curve C : [0, r] — > S 2 . Without loss of 

5 As discussed in the previous section, the choice of a Hamiltonian is not unique, i.e., one can add a 
multiple of the identity operator to Eq. (3.1) without having any physical consequences. 
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generality one can choose a coordinate system in which C does not pass through the south 
pole. This allows one to work with a single coordinate patch of S 2 which excludes the south 
pole. 

In this patch one has pf: 



1 3 

E n (0,0) = bhn, with n = 0, ±- , ±1 , ±- 



E n (9,<p) 

\n;(6,(p)) = e-x j3 e -f J2 e x j3 | n; (0,0)) , 9 e [0, vr), <p e [0, 2vr) 
By definition |n; (0, 0)) are the eigenvectors of H{9 = 0, ip = 0) = 6J 3 , i.e., 

J 3 \n; (0,0)) = ftn|n; (0,0)) . 
In order to compute the operators U(°'(t) and H'(r), one first calculates 

A mn (t) = 4 mn) e\t) + A^tp(t) 



(49) 
(50) 



(mn) 



A 



\(mn) 



d 

(m; (9,^)\ — \n; (8,<p)) 



n 



sm<^ (Ji n )o ~ cos<p (J^ n ) c 
d 



(m; (6,tp)\ — \n; (9,<p)) , 
m(l - cos 0)<5 mn + ^sin0(cos</5 (J^Jo + sin^ (J„ n )o) 



(51) 



(52) 



(53) 



where (J& n )o := (m; (0, 0)| J^|n; (0,0)), (9(t),(p(t)) = C(t) and dot stands for d/dt. In the 
derivation of Eqs. fl52"l ) and (|53|) use is made of the following identities: 

— sin (f J 1 + cos (/? J 2 , 
cos 9 J 3 - sin J 1 , 
cos <£> J 1 + sin <p J 2 . 



e ~ J J z e^ J 



e k " J e « 

-2£j3 7 1 ^2i£j3 



Next step is the calculation of the phase angles a n (t) of Eq. (|7|). These are given by: 

a n (t) = 5 n (t) + ln (t) , 5 n (t) = -btn , ln {t) = -rvy{t) , (54) 

where 



rt ptp(t) 
7(t) := / G?i' (1 — cos 9) cp — / (l-cos0)dy>. 

JO JO 



(55) 
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In (|55"D, ip(0) is set to zero and in the second integral ip is used to parameterize the curve C, 
i.e., 9 = 0((p). 

Eqs. OS), (0), and (g) together with Eqs. @ and @ yield: 



where 



C/(°)(r) = ^e- l[br+7(r)]ri |n;r)(n;0| , 

n 

= J] [#i mn (*) (t)o + ^ m "W (4)o] |m; 0)(n; 0| 

ra/n 



(56) 
(57) 



^ m "(t) := 
flj» n (i) := e 



J{bt+~/](m—n) j 
i[W+7](m— n) / 



sin cp 9 + sin # cos <^ 0) , 
— cos (f 9 + sin # sin </? 0) . 



(58) 
(59) 



In order to further ease the computation of H'(r), one can make the additional assump- 
tion that the energy eigenstates have definite total angular momentum, i.e., 



\J\*\n,(6,<p))=j(j + l) \n,(6,(p)) , n = -j, -j + !,•■■, j . 



(60) 



This assumption is too restrictive for the applications in molecular physics, where one en- 
counters systems with cylinderical symmetry rather than spherical symmetry 0. For a 
magnetic dipole in a classical environment, however, one can safely make this assumption. 



In this case, one can use the well-known relations 11 



J+^J^iJ 2 , J^n; (0,0)) = hC± n \n±l; (0,0)) , C m := - m){j - m + 1) , (61) 
to compute (J^ n )o and (</^,„}o- This leads to 



Substituting Eqs. (^) and carrying out the necessary algebra one has: 



(mn) 



A 



H\t) 



h 

2 



m(l - cos 9)5 mn + - sin6> (e l</, C m £ m „-i + e~ lv C n 5 m _ ln ) 
h\t)Y,C n \n + l;0){n;0\ + h'*(t)J2 C n |n-l;0)(n;0| 



(62) 

(63) 
(64) 
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n 

2 



h\t) e-^ j2 ^ C n |n + 1; (0, 0)> (n; (0, 0) | e<^ + 



%J 2 



/»'•(*) e"^ J £CJn - 1; (0,0)) (n; (0,0) | e^ J 

n 



where 



Re[h'(t)} (cos O J 1 - sin O J 3 ) - lm[h'(t)] J 2 



h'{t) := e l[bt+ ^ ] (i9 + sm9ip) 



(65) 



(66) 



6q := 0(0), and "Re" and "Im" stand for the real part and (— i times) the imaginary part of 
the argument. Furthermore, in the derivation of (|6"5]) use is made of Eqs. (|57D, (|5Ti|), (|(H|). 
V?(0) = 0, and the identities: 



e r J » J = e 1 R J J e ft J ± z J" 1 = cos 9 J 1 - sin 9 J 6 ± z . 



(67) 



Next consider the case in which the magnetic field performs a simultaneous precession 
and nutation, i.e., 

C : (p — tot , cos 9 = cos 9 + i] sm(lujt) . (68) 

Here u is the frequency of precession, and \l\u and rj are the frequency and amplitude of the 
nutation, respectively. Clearly < rj < 1 ± cos 9q and / is an integer. 
The computation of 7 is then straightforward: 



77 

7(i) = (1 — cos 9o)ut — y[l — cos(L;t)] 



Hence for r = T := 2tt/u, one recovers Berry's result: 

7„(T) = -2tt7i(1 - cos0 o ) = -nfi(C) 



(69) 



(70) 



where O(C) is the solid angle subtended by the curve C. A simple consequence of (|70|) is 
that the adiabatic approximation does not detect the effects of the nutation. In fact, one 
has: 



C/(°)(T) = ^e- i[f,T+2 " (1 - coseo)]n |n;0)(n;0| = e -4 J V£ [bT+2 " (1 - cos9o)1J V^ j2 . (71) 
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The first non-adiabatic correction U^(T) to (71) is obtained by integrating H'{t) of 
Eq. fl65|). This involves the evaluation of 

UT) ■= [ T h'(t)dt = f T dt S bt+ ^ ] (i9 + sin 9 0) . (72) 
Jo Jo 

In order to give a closed expression for this integral and subsequently U^'(T), one can 
expand the integrand in powers of 77. This leads to a lengthy calculation which results in: 



Re[/(T)] 



C |sin0 o sin(27r/C) + 
2cos# sin(2vr/C)C 



1 — cos(27r/C) n , 1 „ j. . „ >9x 

-)(-^r + cot0oC-sm^ o C 2 ) 



l-(/C) 2 sin^o 



Im[/(T)] 



sin 3 9 [l-4(l() 2 

C j sin O [1 -cos(2tt/C)] - 

2cos# [l - cos(2tt/C)]C 
sin 3 O [1 -4(/C) 2 ] 



(Zr?) 2 f + 0(^ 2 ) 
1 



sin #0 



-cot0 o C + sin0oC)( 



2 sin(27r/C, 



l-(ZC) 2 ' 



(/^) 2 + 0(/r? 2 ) 



lr]+ 
(73) 

lrj+ 
(74) 



where ( := to / (b — cu cos#o)- I n terms of Re[J(T)] and Im[/(T)], U^'(T) is expressed in the 
form: 



f/ {1) (T) = -i £/ (0) (T) {Re[/(T)] (cos 9 J 1 - sin 9 J 3 ) - Im[/(T)] J 2 } 



.1 e ^J\-iW+2, { i-cose o) ].fi { Re [/(T)] J 1 - Im[/(T)] J 2 } e^ j2 . (75) 



In particular, consider the special case, 9q = tt/2. Then £ = and Eqs. (fHD and 
reduce to: 



[/(o)( T ) 
Re[J(T)] 

Im[/(T)] 



C sin(27r/C) + 



1-C 2 

■1-00 



r )[l-cos(27r/C)] 



CI- cos(2tt/C) + 



(-7 1 T7^)sin(2vr/C) 



i-(/C) 2 ' 



/17J + 0(/r/ 2 ) , 
fyl + G(/r/ 2 ) 



< -•■'( ■ -■' •( - {iWi'/'j ./ ; hn 7i7 ); ./ 2 } ( '- 



ft 



The adiabaticity parameter v of Eq. (|T2"D can also be easily calculated in this case: 

v = cJl + {lr j YC + 0{lr j 2 ) , 



(76) 
(77) 

(78) 
(79) 

(80) 
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where c := C-j/2 = yj(j + 1/2). The appearance of ( in this expression, is an indication 
of the fact that the leading order term in U^(T) is of order v. In fact, for the two extreme 
cases: lr\ <C 1 and Irj ^> 1 (with Irf <C 1 so that one can neglect 0(lrj 2 )) this is manifestly 
seen. 



As indicated by Eqs. (|76]-[FjJ), the nutation of the direction of the magnetic field only 



contributes to the non-adiabatic part of the evolution operator. 



7 Summary and Conclusion 

In this article, I have tried to address some of the basic issues regarding the meaning of the 
adiabatic approximation and the implications of their resolution to the phenomenon of the 
geometric phases. A precise definition of an adiabaticity parameter, although being rather 
unimportant in practice, provides an objective criterion for the applicability of the adiabatic 
approximation. It also allows one to deal with some of the misconceptions concerning the 
subject. Specifically, it is shown that there is no relation between the length of the duration 
of the evolution and the validity of the adiabatic approximation. Furthermore for periodic 
Hamiltonians, the relation between the adiabaticity of the evolution and the period of the 
Hamiltonian is shown to be not a simple proportionality. For example for the system de- 
scribed by Eqs. ( |48| ) and (|68]), both the period T and the nutation parameter lr\ might be 
arbitrarily large, in which case the adiabatic approximation is not valid. An example of 
a periodic system with an arbitrary small period for which the adiabatic approximation is 
exact, is a system with constant energy eigenstates but time dependent energy eigenvalues. 
A typical example of such a system is a magnetic dipole in a magnetic field whose magnitude 
is rapidly changing but its direction is fixed. Clearly, in this case, the adiabaticity parameter 
v vanishes identically and the adiabatic approximation is valid exactly. 

A perturbative expansion of the evolution operator U (r) is also proposed. This expansion 
generalizes the adiabatic approximation and computes the non-adiabatic corrections to the 
adiabatically approximate expression for U(t). Furthermore, it separates the adiabatic and 
non-adiabatic parts of U(t) and yields the non-adiabatic part as a time-ordered exponential 
defined by a Hermitian operator H'(t). In a basis which diagonalizes the initial Hamiltonian, 
H'(t) is off-diagonal. Moreover, the dominant contribution to this operator and consequently 



17 



to the non-adiabatic part of the evolution operator comes from the nearest energy levels. 
Particularly interesting is to view H'(t) as a transformed Hamiltonian. This corresponds 
to a time-dependent unitary transformation of the state vectors which undoes the adiabatic 
part of the evolution operator. 

The subject of the general time-dependent unitary transformations of the state vectors 
has also been reviewed and a simple application of such a transformation is used to yield 
a splitting of the total phase of a cyclic state into its dynamical and geometric parts for a 
periodically changing Hamiltonian. This is followed by a discussion of physically equivalent 
quantum systems and the conclusion that unlike the dynamical (total) phase of a cyclic state, 
the difference of the dynamical (total) phases of two cyclic states is a physical quantity. 

Finally, the theoretical developments have been applied to the analysis of a quantum 
system consisting of a magnetic dipole subject to a magnetic field. In particular, an explicit 
formula for the operator H'(t) has been derived for the general case and the case of a 
precessing and nutating magnetic field has been studied in detail. 
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